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1. Introduction 

Experimentally and theoretically, the cleanest environment to measure parton distribution 
functions (PDFs) describing the partonic structure of nucleons is the Deeply Inelastic Scattering 
(DIS) process where the target hadron is probed with a lepton beam. There has been a long history 
of such experiments starting from fixed target electron-proton collisions at SLAC culminating to the 
HERA collider which collided a 27.5 GeV electron/positron beam with a 920 GeV proton beam. 
For the Large Hadron Electron Collider (LHeC) [1], the plan is to use the proton/ion beams from the 
LHC and collide it with an electron beam from a new accelerator providing 50 — 100 GeV electrons. 
In the future the new electron beam could also be used in conjunction with the hadron beams from 
the projected Future Circular Collider (FCC) to further increase the collision energy. 

The DIS measurements at HERA provide precise constraints for the proton PDF analyses 
down to X ~ 10^^, where x is the momentum fraction of the parton with respect to momentum 
of the hadron. However, for nuclear PDFs (nPDFs) there are no collider DIS data available. As 
the kinematic reach of the fixed target experiments is limited to x > 0.005, the current nPDF fits 
[2, 3, 4, 5] are roughly unconstrained below this x-value. This means that the baseline for heavy-ion 
physics at the LHC is not under good control at the moment. The data from proton-lead collisions at 
the LHC provide some further constraints for the nPDFs but even these new data do not provide 
sensitivity to the small-x region. Here we will show, by performing a new global analysis using 
pseudodata with realistic uncertainties, how the LHeC measurements would improve the precision 
of the nPDF analyses. 

In addition to the nPDF studies, the small-x region is interesting also because a breakdown of 
the linear DGLAP evolution is expected when high-enough parton density is reached. This and 
other physics opportunities at the LHeC in electron-ion mode are briefly discussed as well. 


2. Nuclear PDFs 


The DIS cross sections can be written in terms of structure functions Fj{x,Q^) (j = 1,2,3) 
describing the structure of the target hadrons. In the collinear factorization approach the observed 
modifications of the structure functions with nuclear target are absorbed into process independent 
nPDFs ff (x, Q^). These nPDFs can be extracted from data by performing a global QCD analysis. 
There are two different ways to do this: one can either parametrize the ff^{x,Q^) directly (e.g. 
nCTEQ [3]) or one can take some proton PDF set as a baseline and parametrize the nuclear 
modification R^(x, 2^) (e.g. EPS09 [2]) defined as 






( 2 . 1 ) 


where fi{x,Q^) is the free proton PDF.^ In both cases the parameters are fixed by evolving the 
parametrized nPDFs by the DGLAP equations from an initial scale Qq and comparing the calculated 
cross sections to the available data. The accuracy of the globally analyzed nPDFs are thus dictated 
by the precision and kinematic reach of the included data. 

* A set of proton PDFs is required also in the former option in order to compute the nuclear modifications in Fj{x, Q^). 
Thus, there is no fundamental difference between the two. 
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The current nPDF fits are mainly constrained by fixed-target DIS and Drell-Yan dilepton 
production data. Additional constraints are obtained from inclusive pion production in d-t-Au 
collisions at RHIC (EPS09, nCTEQ and DSSZ [4]) and neutrino DIS (DSSZ). Eigure 1 shows the 
kinematic reach of the data in (x, 2^)-plane used in the EPS09 analysis^ and figure 2 shows the 
expected coverage of the DIS measurement at the LHeC. The increase of the kinematic range is 
more than three orders of magnitude both in x and which would be a huge improvement to the 
prevailing situation. 



Figure 1: The kinematic coverage of the different 
data sets nsed in EPS09 global analysis. The dashed 
line shows the applied lower cut in ^ in the analysis. 
Figure from Ref. [2]. 


Figure 2: The expected kinematic coverage of the 
LHeC nuclear DIS data compared to existing fixed- 
target nuclear DIS data. Figure from Ref. [1]. 


The potential of the LHeC DIS data can be quantified by generating a set of pseudodata and 
performing a re-analysis using this set on top of the existing data, and comparing the results with 
the previous fits. The pseudodata are generated according to the theoretical expectations including 
realistic estimate for experimental uncertainties. Here samples of neutral current DIS reduced 
cross section areduced were generated within 10^^ < x < I and 2 < 2^ < lO^GeV^, for more details 
see Ref. [6]. The ratios of areduced between the lead and proton targets are shown in figure 3 for 
different values up to 500GeV^ as a function of x, both for the generated pseudodata and for the 
theoretical expectation from a baseline fit (an EPS09-style fit) including the uncertainties. Eigure 4 
shows the same data but now compared to the theoretical calculations after including the pseudodata 
into the fit. The inclusion of this pseudodata set clearly improves the precision of the theoretical 
result. 

To see the effect of LHeC pseudodata on the nPDE uncertainties we can compare these 
uncertainties before and after the inclusion of this dataset. This comparison is shown in figure 
5 for valence quarks, sea quarks and gluons at the parametrization scale 2o = 1.69GeV^. Due 
to the lack of constraints at small values of x, especially the gluon uncertainties are large in this 

^BRAHMS data have not been included in any analysis. 
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Figure 3: The ratio between the reduced neutral current DIS cross section with lead and proton target using 
60GeV electron beam. The generated pseudodata are shown in red together with the estimated uncertainties 
and baseline-fit predictions in black. The blue bands show the baseline fit uncertainty. Figure from Ref. [6]. 
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Figure 4: Same as figure 3 but after including the pseudodata into the ht. Figure from Ref. [6]. 


region in the baseline fit. However, when the LHeC pseudodata are included these uncertainties are 
significantly reduced. Also the precision of the sea quark nuclear modification is greatly improved at 
X < 0.01. For the valence quarks the effect is not that significant as the large v-region is already well 
constrained and as the quark distributions at small-v are dominated by the sea quark contribution. 
However, one should keep in mind that the fit funcions in these fits are rather restricted (those of 
EPS09) and do not allow much freedom in the small-v region. In principle the current data would 
allow a very different small-v behaviour but, as there are theoretical arguments why one expects to 
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Figure 5: The nuclear modification of the valence quarks (left), sea quarks (middle) and gluons (right) and 
their uncertainties from the baseline fit (blue band) and from the fit including the LHeC pseudodata (red bars). 
Figure from Ref. [6]. 


have the shadowing effect, it was built into the EPS09 fit functions (used in both fits discussed here). 
When there will be data in this region one should allow an adequate amount of freedom in the fit 
function to reduce the bias caused by too rigid fit functions. This extra freedom affects the shape of 
the small-x behaviour and can also increase the uncertainties. For preliminary studies about this 
topic see Ref. [7]. 

Another future improvement would be the inclusion of charged current DIS that would help 
to constrain the flavor dependence of nuclear modifications. Currently there no data that would be 
sensitive to flavor separation^ so usually the nuclear modifications are taken to be flavor independent. 
The new data from production in proton-lead collisions from CMS [8] may provide some handle 
for the flavor separation but the LHeC charged current DIS would be able to constrain this more 
decisively [9]. 

3. Other electron-ion physics opportunities 
3.1 Small-x physics at the LHeC 

As the LHeC is designed to be capable of probing x values down to 10^®, it would be an ideal 
tool to study saturation effects. The linear QCD evolution is expected to break down when the 
partonic density reaches the point where individual nucleons start to overlap with each other as 
illustrated in figure 6. This phenomenon is typically referred to as a saturation effect, see a recent 
overview in Ref. [10] for more detailed discussions. The onset of this effect is characterized by 
the saturation scale which is expected to roughly behave d& . Due to the nuclear mass 

number dependence the saturation effects are more pronounced for heavy-ion targets. So far there 
has not been clear evidence for these effects at the LHC but the clean measurement and the broad 
kinematic reach would make electron-ion collisions at the LHeC a very promising environment for 
searching these effects as indicated in figure 2 that shows the expected Qs(x) with a Pb target. 

^The use of neutrino-nucleus DIS data for this purpose is complicated by the facts that already the free proton fits use 
these data and that the data are available only as absolute cross sections (not ratios). 
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Figure 6; A schematic diagram of the behaviour of the expected saturation scale in ^ and x. Figure from 
Ref. [1], 


3.2 Hadronization in nuclear medium 

As hadron production in proton-lead collisions at the LHC has turned out to be more complicated 
than anticipated [11, 12], it would be important to study a more cleaner environment in order 
to benchmark the hadron production mechanism on nuclear targets. For low energy hadrons, 
hadronization can happen inside the nucleus, possibly leading to absorption before the hadron is 
formed. For higher energy hadrons, hadronization is likely to happen outside the target but still 
some amount of the partonic evolution can take place inside the nucleus. An attempt has been made 
to quantify these kind of modifications by performing a global analysis of nuclear fragmentation 
functions in Ref. [13], but more detailed measurements are necessary to confirm these effects. 


4. Conclusions and Outlook 

Electron-ion collisions at the LHeC, where an unprecedented energy will be reached, have 
a huge potential to reduce the baseline uncertainties for heavy-ion physics. Perhaps this is most 
apparent in the case of nuclear PDFs. High-precision electron-ion collider data would bring the 
nPDF analysis to a similar level of accuracy as free proton PDFs are at the moment. It would also 
help to relax some assumptions in the fits, e.g. flavor dependence, that are currently necessary as 
they are not well restricted by the data. One should also study how the form of the fit function 
affects the uncertainty estimates to be able to quantify the realistic improvements that the LHeC 
would provide. In addition to the nPDFs, there are also other phenomena that are currently under 
discussion such as saturation at small values of x and hadronization in nuclear environment. Both 
could be accurately studied with electron-ion collisions at the LHeC. 
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